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A general expression is given for the change in free energy when a charge tunnels through a 
junction in a one- dimensional array of TV metallic islands with arbitrary capacitances and arbitrary 
background charges. This is used to obtain expressions for the (average) threshold voltage of the 
Coulomb blockade for a few characteristic geometries. We find that including random background 
charges has a large effect on the TV-dependence of the threshold voltage: In an array with identical 
junction capacitances C and gate capacitances C g , the threshold voltage, averaged over the back- 
ground charge, is proportional to TV a , where a crosses over from | to 1 when TV becomes larger than 

2.5y/ C/C g . 



PACS numbers: 73.23.Hk, 73.40. Gk 



I. INTRODUCTION 

Since the pioneering work by Gorter in 1951 single 
charge tunneling effects have -been extensively studied 
in various kinds of geometriesH Research on single elec- 
tronics has-Jed to potential applications in e.g. current 
standardsjju piiltradense integrated digital electronics 
thermometry,™?! and room-temperature memory.13 In 
many of these applications, tunneling occurs through a 
large number of junctions in series. Mps,t theoretical work 
has assumed homogeneous arrays. QY3 The problem is 
that the number of available states at a finite current 
rapidly increases with the circuit size, so that one cither 
restricts the analysis to, homogeneous arrays or adopts 
a numerical approach.li-3 Using modern techniques, it is 
possible to fabricate arrays of metallic islands separated 
by tunnel junctions with almost uniform capacitances. 
It is however very difficult to avoid non-uniform back- 
ground charges on the islands. This is relevant, since 
the charging energy is very sensitive to the background 
charge. 

The aim of this article is to provide results for inhomo- 
geneous one-dimensional arrays of metallic islands. The 
inhomogeneity can be both in the junction capacitances 
and in the background charges on the islands in the array. 
In particular, we study the threshold voltage for charge 
transport. The results obtained are exact within the. 
classical (orthodox) model of single electron tunneling,t3 
which is accurate when quantum size effects and macro- 
scopic quantum tunneling effects may be ignored. 

Using a general expression for the inverse capacitance 
matrix, we calculate in Section || the change in the free 
energy of an TV-junction array due to an arbitrary tun- 
neling event. In Section [II, we focus on the threshold 



voltage for transport Vt, which is an observable quan- 
tity. We find that inhomogeneity of the junction ca- 
pacitances C has a small effect on the threshold voltage 
in large arrays: The expectation value as TV — > oo for 
the threshold voltage of an array without gate coupling 
(gate capacitance C g = for each junction) and without 
background charges is (Vt) — |TVe(C _1 ), with (C _1 ) 
being typically not much different from 1 / (C) . How- 
ever, as we show in Section IV, a random variation in 



background charges may change the threshold voltage 
considerably: In a short array with weak gate coupling 
(N 2 C g /6.25C < 1) and random charges on all TV islands, 
we find (Vt) oc \/~N. In a long array with strong gate 
coupling (TV 2 C 9 /6.25C > 1, but still C g < C), we find 
(Vt) oc TV. We compare our results with experiments. E3 



II. FREE ENERGY 
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FIG. [j]. Schematic diagram of a one-dimensional array of 
TV tunnel junctions. Island i is coupled to island i + 1 by a 
tunnel barrier with capacitance Cj+i, and to a gate electrode 
by an insulating barrier with capacitance C s> j. The capaci- 
tance Ci (Cjv) denotes the coupling of the first (last) island 
to the emitter (collector) electrode. 
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The system under consideration is shown schematically 
in Fig. [|. Within the orthodox model, the state of the 
system is described by the numbers n, of electrons on 
the i-th island, which wc combine in a vector: n = 
(m, i%2, ■ ■ ■ , njv— i)- The tunneling rate, Tk{n), corre- 
sponding to a single e^ctron tunneling from island k — 1 
to island k is given by 



r fc (n) 



AG fe (n) 



e 2 i? fe [l-exp(-AG fc (n)/fc B r)]' 



(1) 



Here Rk is the resistance of the fc-th tunnel junction and 
AGfc (n) is defined as the difference in free energy of the 
final and initial states. The free energy comprises the 
electrostatic energies of the charged capacitors in the sys- 
tem, as well as the potential energies of all electrodesH 
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(2) 



the electrochemical potential of island 
V c ), and by Q e , Q c , and Q g ,i the 



We denote by 4 
i [<j)Q = V e and 
charges on the emitter, collector, and gates, respectively: 



>N 



Q e = Ci(V e ~ fa) + en e , 
Q c = C N (V C - 4>n-i) + en c , 

Qq.i — C gi i(Vg^i — 4>i). 



(3a) 
(3b) 
(3c) 



Here n e (n c ) is the number of electrons that has tun- 
neled from the emitter (collector) electrode through the 
first (last) capacitor. 

The difficulty in determining the energy difference 
AGfc(n) lies in the determination of the electrochemical 
potentials <j) = (<f>i, 02, ■ • ■ , 4>N-i)- They follow from the 
condition that the total capacitive charge on each island 
i equals erii plus a background charge Qqj: 



C g ,i{4>i - V g s 



Ci((pi - <f>i-i) + C i+ i{ 
erii + Qoi, i= 1,2, ...,N— 1. 



(4) 



The background charge Qo,i € (— e/2, e/2) is due to in- 
completely screened charges in the environment of the is- 
land. Eq. ([|) can be written in matrix form as C(f> = Q' . 
with 

Cij = 5ij(Ci + Ci+i + Cg^) — Si+ijCj + Sij + iCi 7 (5a) 
Q'i = erii + Q ,i + CgjVg^ + 5i_iCiV e + Si_ N -iC N V c . 

(5b) 

The capacitance matrix C can be inverted exactly. The 
elements Rij of the inverse capacitance matrix R = C -1 
are given by 



Ri,j — Ci+iCi+2 ■ ■ ■ CjDi^iDj + iD N 1 _ 1 , i < j, 



R 



R 



-1,3- 



(6) 



Here we have introduced the subdeterminants Di (D^-i) 
of the upper left (lower right) capacitance submatrix of 
dimension i. These can be found recursively from 

Di = (C, + d+i + C S)i )A-i - Cf A-2, (7a) 
Di = (C< + C i+1 + C g ,i)D i+1 - Cf +1 D i+2 , (7b) 
D = D N = 1. (7c) 

For a homogeneous array with identical capacitances, 
C\ = Ci = . . . = Cn and C g> \ = C g ^ = ■ ■ ■ = C„.n-i, 
we recover the inverse capacitance matrix of Ref. |l2] . 

We now derive a general expression for the difference 
in free energy AGfc(n) when an electron tunnels from is- 
land k — 1 to island k. Applying Eq. (Q) and making use 
of the orthogonality relation 

(G, + Gj+x + C g s)Rij = CiRi^ij + Cj+i-Rj+xj + Sij, 

(8) 

we find that AG&(n) takes the form 



J 



N-l 



AGfc(n) = — (Rk-l,k-l + Rk,k — Rk-l,k — Rk,k-l) + e Qi(Ri,k-l ~ Ri,k) 



N-l 



+ e(V e - V g ,i)A hk + ( v 9,i-i - Vg,i) A i,k + e{V g ,N-x - V c )A N , k , 

i=2 

Ai,k = Ci(Ri-i.k + Ri,k-i — Ri-i,k-i — Ri,k) + SiM- 



(9a) 
(9b) 



Here, Ri,n = Ro,i = is implied. 

Although we are now able to construct all relevant 
transition rates from expressions (^|) and (^|), the ana- 
lytic evaluation of the current-voltage characteristic at 
arbitrary voltage remains a technically involved problem. 



The threshold voltage, however, is determined by a single 
transition rate and is therefore easier to evaluate. In the 
next two sections, we apply our results to this quantity 
for several characteristic geometries. 
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III. THRESHOLD VOLTAGE 

Electron transport through a one-dimensional array is 
realized by a sequence of tunneling events through all 
junctions between the emitter and the collector (we refer 
to this as a tunneling sequence). At zero temperature, 
a specific tunneling sequence contributes to the conduc- 
tance if the free energy difference of each tunneling event 
in the sequence is positive. The threshold voltage V t of 
the Coulomb blockade is the smallest voltage at which a 
current can flow through the array at zero temperature. 
When \V e — V c \ < \V t \, there exists no conductive tunnel- 
ing sequence. We first consider the simple case where the 
system is not gated (C g ^ — for all i), and then discuss 
the turnstile configuration, i.e. an array which is coupled 
to a gate electrode via a single island: C g ^ = C g Si^ n . 



The threshold voltage is determined as follows. For an 
initial charge state, we determine the minimal activation 
energy eVt.k(q) to allow a tunneling event in the array, 
as well as the corresponding final charge state. The final 
charge state becomes the initial state in the next step. 
The minimal activation energy for the new charge state 
and the corresponding final charge state are again de- 
termined, and this procedure is repeated until one elec- 
tron has been transported from emitter to collector. The 
largest of the activation energies found equals eVt- In the 
special case that all background charges are zero, one has 

V t = |e (j^ 1/d - Max[l/d, 1/C 2 , . . . , 1/C N ]J , 

(13) 



A. No gate coupling 

In the absence of gate coupling, the determinants D 
and D, following from Eq. (]?]), have a simple form. For 
convenience, wc introduce the notation 



which is an extension of the cesult Vt = 
±eMin[l/Ci, 1/C 2 ] for a double junctionJl3 For N -► oo, 
Vt has a Gaussian distribution with average iA^e(C _1 ) 
and variance VarVt = -|iVe 2 VarC _1 . 



i=k+l 



S'' 



In terms of these quantities, 

D k = C 1 C 2 ---C k+1 S k+1 , 
Dk = CkCk+i ■ ■ ■ CatS^-i, 
Rij = S i S j /S N ,i<j. 



(10) 



(11a) 
(lib) 
(11c) 



We further define q 
e ~ (Qo,i,Qo,2, ■ ■ ■ ,Qo,n-i) 



n + q , q = 
From the condition 
AGk(q) = 0, we determine the threshold voltage V ti k(q) 
for tunneling through capacitance C k at arbitrary occu- 
pation q of the array: 



k-l 



AT-l 
i—k 



QiSi. (12) 



J 



B. Turnstile configuration 



We next consider a turnstile configuration, i.e. an ar- 
ray with a single gate electrode coupled capacitively (ca- 
pacitance C g ) to island n. The elements of the inverse 
capacitance matrix are then given by 



Ri 
Ri 



(S l + C g S n S l n )S J (S N + C g S n S n )-\ n<i<j 



N 



Con q 



i < n < j 



Ri,j S {Sj + CgSj S n )(S + CgS S n ) 



Rj,i — Ri,j • 



i < j < n 
(14) 



In order to determine the threshold voltage Vt t k{o), we 
have to distinguish between k < n and k > n. From Eqs. 



<® and (h 

voltage 



-we find that Vt,k(q) now depends on the gate 
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k < n 
k > n 



(15) 



where 5" is defined as in Eq. ( |l0|) in terms of modified 
capacitances C: 



Ci(l 
C' l = C l (l 
C' l =C l (l 
C[ = C,(l 



\CgS n )(l 

2~CgS n ), 

^CgS n ), 

\C g S n ){l 



CgS V 



c g s n y 



k < n,l < n, 
k < n,l > n, 
k > n,l < n, 
k > n.l > n. 



(16) 



IV. BACKGROUND CHARGE 



The background charge in a single-electron tunneling 
device has a large influence on its properties. For exam- 
ple, by tuning the background charge in a double junc- 
tion with one gate one can set the threshold voltage to 
any value between zero and e/(2C + C g ). In this section, 
we investigate the effect of background charges on the 
threshold voltage of an array of tunnel junctions. For 
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reasons of clarity, we choose identical junction capaci- 
tances in the following (Cj = C for all i). We start by 
investigating an array with a non-zero background charge 
on a single island. We then give ensemble-averaged re- 
sults for random background charges on all islands and 
compare with the experiments of Dclsing et al£3 

In the absence of gate coupling (C S) , = for all i) 
and for a non-zero background charge go,™ = Qo,m/e on 
island m, there are three initial tunneling events which 
may form the bottleneck for conduction: 

• transfer of an electron from the emitter to the first 
island (electron injection through junction k = 1); 



• tunneling through junction k = 
through junction k — to if go,-, 
creation at island to); 



m+1 if qo. m > or 
i < (electron-hole 



• transfer from the last island to the collector (hole 
injection through junction k — N). 

An analysis of the corresponding tunneling sequences re- 
sults in the threshold voltage 



<Z0,m > 0. 
<?0,m < 0. 



1 - 2Min[rTO7o, m , (N - m)(l - q , m )}) , 

(17a) 

1 - 2Min[m(l - \q , m \), {N - m)\q , m \]) , 

(17b) 



For a uniform distribution of qo. m between ±i and a 
uniform distribution of to between 1 and N — 1 its ex- 
pectation value is (Vf) = (5N — 7)e/12C, with variance 
VaxVt = (e/2C) 2 (A^ + l)(3^ 2 -5A^ + 8)/180^. The ex- 
pectation value is slightly smaller than for a homogeneous 
array without background charges: V t = (N — l)e/2C. 
In the limit N — > oo the root-mean-square deviation is 
rmsVt cx Ne/C, of the same order as the threshold volt- 
age itself. 
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FIG. g. Derivative of the average threshold voltage with 
respect to the array length N, for ensembles of arrays with 
identical capacitances (Ci = C and C 3t i = C g for all i) and 
random background charges on all islands, calculated from 
Eq. (|l8|). The average is determined numerically from ensem- 
bles of 10000 samples for N < 128 and ensembles of 1000 
samples for larger arrays. A cross-over from (Vt) oc N 1 ' 2 to 
(Vt) oc TV occurs at N c » 2.5^/C/Cq. Solid lines are the 
extrapolation formulas (G2|) and (E3h. The dashed curves are 
obtained from the result (hu) for zero background charges and 
Cg/C — (upper curve) and C g /C — 0.01 (lower curve). 

We next consider a one-dimensional array of equally 
gated islands (d — C, C g .i — C g for all i). In Refs. Q 
and |l^ the charge transport in homogeneous arrays by 
soliton-like excitations was introduced. In terms of the 
soliton width A" 1 = [2arsinh v /C 9 /4C]- 1 of Ref. |, the 
threshold voltage for an electron tunneling through junc- 
tion k is given by 



JV-l 



VtA<f) = ^ - 2 I](* + la) sinh(zA) cosh[(/V -k+ |)A] + 2 ^ ( ft + q g ) sinh[(/V - i)X] cosh[(A: - i)A] 

\ i—1 i—k 

+ sinh[(/V - |)A] - cosh[(/V - 2k + 1) A] sinh ^ J (^sinh A cosh ^ cosh ^~ 2 2 fc + 1 ^ A 



(18) 



Here, the gate-induced charge q g = C g [V g — \ {V e + V c )] 
acts as an offset on the background charge. The average 
threshold voltage (averaged over the background charge) 
is therefore independent of V g . For N = 2, we find 



(V t ) = e/(AC + 2C g ). 



(19) 



In the absence of background charges and for q g = 0, 
we find 



V t = 



e sinh[(AT- l)A/2] 
2C cosh(7VA/2) s inh(A/2) : 



(20) 



which approaches a constant value as N — > oo, provided 
A ^ 0, i.e. provided C g /C ^ 0. In Figure | we show 
the effect of random background charges on all islands 
in arrays of different lengths for several gate couplings, 
as calculated from Eq. (|l8|). The averages are computed 
numerically by putting a random charge go.fe G ( — §j |) 
on each island k. The dependence of (Vt) on the ar- 
ray length differs drastically from the result (^0|) without 
background charges: Instead of a threshold voltage which 
exponentially approaches a constant value as N — *■ oo, we 
find (Vt) oc y/N — 1 for small arrays, with a cross-over to 
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a linear TV-dependence for large arrays. For C g <C C, the 
array length N c at which the cross-over occurs is found 
to be 2.5 times the soliton width, 



N r m 2.5 



\fc/c 9 



2.5A" 



(21) 



For N < N c the average threshold voltage is well de- 
scribed by an extrapolation of the result (fL9h for N = 2: 



(Vf) = 



4C + 2C g y/2-l' 
For N > N c we can describe the numerical data by 



(22) 



(V>) = (V<) N=Nc + (N- N c ) 



d(V t < 



1 



iC + 2C g y/2- 1 



dN 
N + N, 



N=N C 



- 1 



(23) 



The cross-over to a linear iV-dependence supports the 
intuitive idea that the background charge in the array 
is screened beyond N c . The rms deviation rmsVt = 
0.31e(v^V - 1)/(2C + C g ) for all N. The rms deviation 
of the threshold voltage for tunneling through a specific 
junction k has a much stronger dependence on N than 
rmsVf itself: rmsV^fc oc N 3 / 2 . Since V t is chosen as the 
maximal threshold voltage in a sequence of N minimal 
values for single tunneling events, the fluctuations in V t 
are smaller than those in Vt k ■ 
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FIG. q. Comparison of experimental threshold voltages 
(taken from Ref. [p| solid dots) with the result of Eq. ([lq), 
averaged over the random background charge (open squares 
with error bars) . We used identical gate and junction capac- 
itances, with Cg/C = 0.044 (N c = 12), as estimated in Ref. 
|l6| . There are no adjustable parameters. 

In Figure [| we compare the threshold voltage from 
Eq. (|l8|), averaged over all background charges, with 
experimental threshold voltages for arrays of different 
lengths.^ We used the values C = 0.28 fF and C g =0.012 
fF from Ref. [L6L giving N c = 12. Thus, the experimental 



results are in the regime of a linear dependence of (14) 
on N. The qualitative agreement is satisfactory, without 
any adjustable parameters. 

In conclusion, we have derived an exact analytical 
expression for the threshold voltage Vt } k(q) for tunnel- 
ing through a junction k in a one-dimensional array of 
N metallic islands at arbitrary occupation q of the is- 
lands. We have calculated the average threshold volt- 
age for transport and its fluctuations in a few simple 
cases. In particular, we have found that including ran- 
dom background charges results in a N a dependence of 
(V t ), with a = i for iV < 2.b^/CjCTg and a = 1 for 
N > 2.5-JC/Cg. We have made a comparison with 
the available experimental data on gated one-dimensional 
arrays,c!j and found a reasonable agreement. 
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